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ABSTRACT 

The delay in light travel time along the line of sight generates an anisotropy in the 
power spectrum of 21cm brightness fluctuations from the epoch of reionization. We 
show that when the fluctuations in the neutral hydrogen fraction become non-linear 
at the later stages of reionization, the light-cone anisotropy becomes of order unity 
on scales > 50 comoving Mpc. During this period the density fluctuations and the 
associated anisotropy generated by peculiar velocities are negligible in comparison. 
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1 INTRODUCTION 


Fluctuations in the 21cm brightness from cosmic hydro¬ 
gen at redsliifts z > 6 were sourced by the primord ial 
densi t y perturbations f rom in flation dLoeb fe Zaldarriaeal 
120041 ; IBarkana fe Loebl |2005a|B) as well as by the radi¬ 


ation fro m galaxies dScott_fc_Rggs| 1990k_ Madau et al.l 


1991]; _ [Furlunetto ^Zaldarriaga. fe Hernauist 20041 ; 

Barkana fe Loebl l2005bfh These two different compo¬ 


nents can be separated based on the angular dependence of 
the 21cm fluctuati on power spectrum which is induced by 
peculi ar velocities llBharadwai fe A li 2004; Bark ana fe Loebl 
I2005af ). Uncertainties in the cosmological parameters lead 
to an additional apparent anisotropy due to the Alcock- 
Paczynski effect; after accounting for constraints from 
the cosmic microwave background, this effect can still 
produce up to a 10% anisotropy which may be_identifiable 
because of its particular angular structure dNusse 1 12005| : 
IBarkana! 12006)') . Several low-frequency arrays that could 
potentially detect the redshifted 21cm signal are being 
built around the globe, including the Primeval Structure 
Telescope (web.phys.cmu.edu/~past), the Mileura Widefield 
Array (web.haystack.mit.edu/arrays/MWA), and the Low 
Frequency Array (http://www.lofar.org). 

In this work we examine an additional source of 
anisotropy in the 21cm power spectrum between the direc¬ 
tions parallel and transverse to the line of sight. Previous 
discussions ignored the delay in light travel time (i.e., the 
“light-cone” constraint) along the line of sight [but see a 
rough estimate of the effect at the beginning of reionization 
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in an appendix in lMcQuinn et ahl (l200d )l. In particular, al¬ 
though the 21cm power spectrum is expected to be measured 
through averaging over three-dimensional volumes of finite 
radial extent, the power spectrum was previously evaluated 
at a fixed time slice of the universe, ignoring the fact that 
a fixed observing time implies a varying emission time as 
a function of distance from the observer. Here we evaluate 
the amplitude of the anisotropy which is sourced by this 
delay. This time-delay effect is not related to real causal 
effects or light-crossing times. Since two points at a given 
separation are in general seen at different redshifts, the cor¬ 
relation function, averaged over all such points, is affected 
by the change with time of the statistics of ionization (i.e., 
the distribution of H II regions and their correlation with 
the underlying density field). 


2 MODEL 

IBarkana fe Loebl (120041 ') showed that the biased large- 
scale fluctuations in the number density of galaxies at 
high redshift imply that the characteristic bubble ra¬ 
dius is determined by correlated groups of galaxies 
rather than by the bub ble sizes of indiv idual galaxies. 
iFurlanetto. Zaldarriaga, fc Hernauistl (120041) developed a 
semi-analytic model based on this realization, and used it to 
predict the distribution of H II bubble sizes around a given 
point (i.e., the one-point distribution of bubbles). Unfortu¬ 
nately, this model cannot be directly applied to estimating 
two-point distributions. Upcoming 21cm experiments will 
only have the sensitivity to probe statistical measures of 
the 21cm fluctuations and will therefore focus on the three- 
dimensional power spectrum (or, equivalently, the two-point 













































2 R. Barkana and A. Loeb 


correlation function). We formulate here an approximate 
model of the 21cm correlation function in order to estimate 
the time-delay anisotropy. 

Suppose that gas with an overdensity <S and a neutral 
fraction x„ lies in some direction at a distance correspond¬ 
ing to 21cm absorption at a redshift z. Then the resulting 
21cm brightn ess te mperature offset relative to the CMB is 
jMadau et al.lfl997l j 

where we have assumed that the IGM has been heated so 
that T s Tomb, and we have substituted the concordance 
values (jSpergel et al .1120031 ) for the cosmological parameters 
flf,, h, and fi m . In order to calculate the two-point correla¬ 
tion function of the 21cm brightness temperature, we must 
construct the correlated distribution of <5 and x n at two dif¬ 
ferent points observed at two possibly different redshifts. We 
first estimate the characteristic bubble radius R c h at each 
redshift, and calculate the mean ionized fraction of the in¬ 
terior of regions of radius R c h as a function of their mean 
densities. We then calculate the correlation function at two 
points by integrating over the joint probability distribution 
of their overdensities. 

We first must estimate the characteristic radius R c h at 
a redshift z, defined so that for a given point P, the average 
number of ionizing sources within P c h typically determines 
whether P is ionized or not. In estimating R c h we must 
balance two competing trends. On the one hand, if a large 
region contains enough sources to fully reionize itself, then 
it will indeed fully reionize regardless of whether small por¬ 
tions within it contain a low density of galaxies; this simple 
fact favors larger radii over smaller ones. On the other hand, 
the ACDM density power spectrum implies that the typical 
magnitude of density fluctuations decreases with scale, so 
that large regions are unlikely to have a strong positive fluc¬ 
tuation that will increase the number of galaxies sufficiently 
for them to fully reionize. Thus, we estimate R c h as the 
largest radius for which a l-cr density fluctuation results in 
full reionization of a sphere of this radius. 

Towards the end of reio nization, t h is sim pl e est imate of 
P c h approaches infinity, but lWvithe fc Loebl (12004) showed 
that the bubble radius is limited by the finite light travel 
time to ~ 70 comoving Mpc at redshift ~ 6; we generalize 
their argument and set an upper limit on R c h at any stage of 
reionization. In general, in order to have a coherent ionized 
bubble of radius P c h, various portions of the bubble, some 
of which are relatively rich with galaxies and some of which 
have fewer ionizing sources, must effectively exchange pho¬ 
tons between them. There are two important timescales: the 
typical (i.e., l-cr) scatter Scatter in the time when different 
regions of size P c h reach a given stage of galaxy formation (in 
terms of a given number density of galaxies); and the light 
crossing time fn g ht of the radius P c h. Now, if fn g ht > tscatter 
at a given redshift z, then such a large region will not con¬ 
stitute a single, coherent bubble; e.g., a low-density por¬ 
tion of this region that has not produced enough photons to 
reionize itself will indeed not be fully reionized at redshift 
z, since the extra ionizing photons needed from the distant 
portions of the region will typically have been produced a 
time ~ tscatter earlier but requ ire a time ~ fligh t to re ach the 
low-density region. We follow IWvithe fe Loebl d2004 in us¬ 


ing the extended Press-Schechter model dBond et al.l Il99lh 
to estimate this upper limit on P c h. 

The progress of reionization depends on the total col¬ 
lapse fraction of gas in galactic halos, F co i, which is deter¬ 
mined by integrating the halo mass function over all ha- 
los that host galaxies. We use the halo mass function of 
ISheth fe Tormeii dl99S ) whic h fits numerical simulations ac¬ 
curately ( Jenkins et al. 120011 ) . In order to calculate fluctua¬ 
tions in the gas fraction that resides in galaxies within dif¬ 
ferent regions, we adj ust the mean halo distri bution based 
on the prescription of iBarkana fe Loebl |20 q 4 which fits a 
broad range of simulations. As gas falls into a dark matter 
halo, we assume it can fragment into stars only if its virial 
temperature is above 10 4 K, which allows for efficient cool¬ 
ing mediated by atomic transitions. We assume that stars 


dominate over mini-q uasars, as expected at high r edshift 
dWvithe fe Loebl 12003 1 [but see iMadau et aid (12004 ]. The 
stellar spectrum depends on the initial mass function (IMF) 
and the metallicity of the stars. We consid er two exam ples. 
The first is the locally-measured IMF of ISca k] dl998ll with 
a metallicity of 1/20 of the solar value, which we refer to 
as Pop II stars. The second case, labeled as Pop III stars, 
consists entirely of zero metallici ty M > 100 tfp stars, as 
expected for the earliest galaxies (iBromm fe Larsonll2004 . 

If we were to treat each sphere of radius R c h as an 
isolated region, then the mean ionized fraction in a region 
of mean perturbation 5 would simply be related to the gas 
fraction in galaxies F co i: 


Si = ( 1 + F^) Fcol(2,jRch ’^ ’ (2) 

where -Won is the overall efficiency of emission of ionizing 
photons per baryon in galactic halos, the factor of 0.76 con¬ 
verts between hydrogen and baryon number densities, and 
the (minor) addition of 1 counts the gas in galaxies (which 
need not be ionized by the photons that escape from galax¬ 
ies). Note that Won equals the number of ionizing photons 
produced per baryon in stars, multiplied by the star forma¬ 
tion efficiency /, and by the escape fraction / esc of ionizing 
photons from halos. 

In reality different regions are not independent. The 
overdense regions produce more ionizing photons than they 
need to fully reionize themselves, and the extra photons ex¬ 
pand into the surrounding lower-density regions. In order to 
make our model consistent with the total number of ionizing 
photons produced at each redshift, we increase the Xi values 
by the amount necessary to achieve consistency. Now, the 
leakage of additional photons is expected to depend on the 
mean <5 of the target region. In particular, the deepest voids 
are likely to be the farthest distance away from the densest 
regions. We estimate the relative fractions distributed into 
different regions as follows. For each S, we assume that the 
increase in Xi within a region of that 8 is proportional to the 
chance that such a sphere lies next to a sphere with 8 = Si, 
where Si is a high overdensity representative of those regions 
that have overproduced ionizing photons. We adopt for Si 
the minimum value of S needed for a sphere to achieve an 
Xi = 1 on its own, i.e., without the help of surrounding 
regions. In general, we emphasize that the redistribution 
of photons has a modest effect within our model, consis¬ 
tent with our selection of R c h in such a way that photons 
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from larger distances are unlikely to play a major role in the 
statistics of reionization. 

When we calculate the correlated ionization states of 
two points at redsliifts z\ and 22 , we assume that the ion¬ 
ization probability of each point is determined by the pertur¬ 
bations 5 of the surrounding regions of radius R c h{zi) and 
Rch(z 2 ), respectively. However, the 21cm brightness tem¬ 
perature Tb is determined not only by the neutral fraction 
but also by the perturbation <5, which is in general differ¬ 
ent from the spherically-averaged perturbation 5. While 8 
is determined by large-scale power above the scale of R c h, 
the difference <5 — <5 is determined by the additional small- 
scale power within the bubb le. In the spirit of the extended 
Press-Schechter model dBond et al.ll 19911 1 , we assume that 5 
and 5 — 8 are approximately statistically independent. We 
thus add the portion 5 — 5 at each point independently of 
the values of 5 and the neutral fraction a: n of the two re¬ 
gions, with a correlation function for the small-scale portion 
of 8 set equal to the total £(r) of the density held minus the 
smoothed £(r) that determines the joint statistics of 5. 

When calculating the correlated ionization states of two 
points we must carefully interpret the values of x n ,i and x n ,2 
in the regions surrounding the two points as probabilities of 
having neutral gas. For example, the mean value of x 2 at 
point 1 is not a; 2 1; which it would be if the gas were partially 
ionized with a neutral fraction of x n ,i, but is instead x„,i 
since the gas has x„ = 1 with probability x n ,i and x„ = 0 
with probability 1 — x n ,i■ Thus, the mean expected value 
of Xn,iXn ,2 must be calculated as the joint probability that 
both points are neutral. Given the individual probabilities 
x n ,\ and x n , 2 , there is some probability p sam e that the two 
points are ionized within the same bubble. Assuming that 
their ionization probabilities are independent if they are not 
within the same bubble, we find that 

(Xn,lX n ,2) — Xn,lXn : 2/id Psame) - (3) 

Now, in order to estimate p ,« for two points separated by 
a comoving distance r, we note that if one point is ionized 
by a bubble of radius R c h then the probability that this 
(randomly-placed) bubble contains a second point a distance 
r away is zero if r > 2R cbl and 

f(r,R ch ) = l-^(R 2 ch -±r 2 ) (4) 

otherwise dFurlanetto, Zaldarriaga. fe Hernquistl 12004) . 

This yields two estimates of p same , starting from either 
point, so we take the average as our best estimate: 

Psame « ^ [(1 ~ X n ,l ) f(r, Rch (zi )) + (1 - X„, 2 ) f(r, R c h (z 2 ) )] , (5) 

with the constraint that p same must not be higher than the 
separate ionization probabilities of each of the two points. 

The two-point correlation function £ of the 21cm bright¬ 
ness temperature Tb is a function of the distance r between 
the two points and their two redshifts. In order to focus on 
the anisotropy as a function of the angle 6 between the line 
of sight and the vector connecting the two points, we pa¬ 
rameterize £ as a function of r, p = cosd, and the redshift z 
at the midpoint (in terms of comoving distance) of the two 
points. We calculate the expectation value 

£ = ((T b}1 - f b (z!)) x (T b}2 - f b (z 2 ))) , ( 6 ) 

assuming that the mean T b at the appropriate redshift has 


been subtracted from T b at each point. Such a subtraction 
is expected in future observations since the radio measure¬ 
ments take the Fourier transform of the intensity in the x 
and y directions, removing any sensitivity to a varying T b (z) 
as long as the line k x = k y = 0 in Fourier space is avoided. 
We also assume in our results below that T b (zi) has been 
multiplied by [(1 + z)/{ 1 + 21 )] 1 / 2 , and similarly for T b (z 2 ), 
in order to avoid any anisotropy due to the simple redshift 
scaling in equation ©. 

We calculate the correlation function explicitly by aver¬ 
aging over the joint normal distribution of the perturbations 
5 of the regions of radii -R c h(zi) and R c h{z 2 ) surrounding the 
two points. We also calculate the mean value T b (z) at each 
redshift by averaging the quantity [(1 + (5)*n] over the prob¬ 
ability distribution. Note that we do not include peculiar 
velocity fluctuations here since we focus on the late stages 
of reionization and show that in this regime the 21cm fluc¬ 
tuations and their anisotropy are dominated by ionization 
fluctuations. 


3 RESULTS 

We illustrate our predictions for two cases which both as¬ 
sume atomic cooling, a star formation efficiency /* = 2%, 
and an escape fraction of ionizing photons from halos of 
/esc = 10%; however, our low-* case assumes Pop II stars 
(§ 2) and achieves complete reionization at redshift 7.35, 
while our high -2 case assumes Pop III stars and completes 
reionization already at redshift 13.7. Note that the parame¬ 
ter values given here assume that only one ionizing photon 
per hydrogen atom is needed in order to achieve reioniza¬ 
tion of the IGM. If we assumed instead that several photons 
are needed due to recombinations, then this would require 
an increase in the value of /* by the same factor as the re¬ 
quired number of photons per atom. This approximation is 
fairly accurate as long as the number of times each hydrogen 
atom must be reionized is roughly the same in different re¬ 
gions. Indeed, the recombination rate is not expected to vary 
greatly in different large-scale regions, since the fluctuation 
in the mean density of such regions is small. Recombinations 
can be included in a future refinement of our model. 

Figure |T] shows the two reionization histories, as well 
as the characteristic radius R c h and the standard deviation 
of the mean density 8 in a sphere of radius R c t. At each 
redshift the probability distribution of 5 in the spherical re¬ 
gion is a Gaussian with this standard deviation; we have 
not included non-linear corrections which are particularly 
important early in reionization when R c h is relatively small. 
The characteristic bubble radius increases during the epoch 
of reionization as larger groups of galaxies are able to fully 
reionize the region containing them. At the later stages of 
reionization, when the cosmic mean x„ < 10%, R c h is deter¬ 
mined by the light travel time of ionizing photons. 

In Figure [2] we illustrate the correlation function £ of 
21cm brightness temperature late in reionization in the high- 
2 case. If we considered density fluctuations alone (along 
with a uniform x n {z) equal to its cosmic mean value at each 
redshift), then ^(r) would roughly equal the 1 -a density fluc¬ 
tuation on scale r, times the cosmic mean T b , all squared. 
The growth factor and x„{z) both evolve slowly at this stage 
of reionization and produce a negligible anisotropy. Adding 
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Figure 1. Reionization history. We show the cosmic mean 21cm 
brightness temperature (top panel), redshift z (next lower panel), 
the comoving characteristic radius R c |, (next lower panel), and 
the 1-a density fluctuation in spheres of radius R c ^ at z (bottom 
panel), each vs. the cosmic mean neutral fraction x n . Each panel 
shows our low-z case (solid curve) and our high-z case (dashed 
curve). 


in peculiar velocities generates a significant anisotropy, but 
only causes fluctuations that are of order the density fluc¬ 
tuations. However, the ionization fluctuations are far larger 
and indeed non-linear, since while the bubble radius R c h is 
large and the corresponding density fluctuations are small 
(see Figure [T]). various regions can have values of x„ cover¬ 
ing almost the entire possible range of 0-1. In particular, we 
find that although the l-a density fluctuation on the scale 
of R c h ~ 20 Mpc is ~ 0.05 at this redshift, a negative 2-a 
void on this scale has an x„ ~ 60% (compared to a cosmic 
mean x n of 25%), while a positive l-a overdense region is 
almost completely reionized. 

Figure 03 illustrates the 21cm correlation function near 
the end of reionization in the low-z case. The ionization 
fluctuations again dominate, since the deep voids still have 
a substantial neutral fraction that is much higher than the 
cosmic mean x n - The strong anisotropy here is related to the 
rapid reionization of the deep voids at the end of reioniza¬ 
tion; in particular, £ drops to zero when the lower-redshift 
point approaches the redshift where cosmic reionization is 
completed. Although we have predicted the general trends, 
we emphasize that the results at the end of reionization de¬ 
pend on how rapidly the deep voids get reionized by photons 
from denser regions, and thus precise quantitative predic¬ 
tions must await fully self-consistent calculations of large- 
scale radiative transfer. 

Figure [4] summarizes the evolution of the 21cm correla¬ 
tion function during the latter half of the reionization epoch. 
Early in reionization, ionization affects only the rare high- 
density regions, and the anticorrelation of x„ and 5 leads 
to a low £ even though the cosmic mean brightness tem- 



Figure 2. Correlation function of 21cm brightness temperature 
as a function of r or //, at a late stage of reionization (x n = 0.25) 
by Pop III stars. The left panel considers comoving r = 20 Mpc 
(solid curve), 50 Mpc (short-dashed curve), and 100 Mpc (long- 
dashed curve). The right panel considers angle /r = 0 (solid curve) 
and fi = 1 (dotted curve). Actually shown is |£|, since £ becomes 
negative when r ~ 200 Mpc and fi = I. Also shown are the predic¬ 
tions of a model with density and peculiar velocity fluctuations 
but no ionization fluctuations, for fi = 0 (short-dashed curves) 
and fi = 1 (long-dashed curves). 



Figure 3. Correlation function of 21cm brightness temperature, 
same as Figure [2] except near the end of reionization (x n = 0.05) 
by Pop II stars 
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Figure 4. Correlation function of 21cm brightness temperature 
as a function of the cosmic mean neutral fraction x n . Each panel 
corresponds to a comoving distance r as indicated. We consider 
Pop II stars and an angle of p, = 0 (solid curves) or 1 (dotted 
curves), or Pop III stars and an angle of p = 0 (short-dashed 
curves) or 1 (long-dashed curves). 


perature is relatively high. Later on, as reionization spreads 
to more typical regions, <5 fluctuations produce amplified x n 
fluctuations and £ increases, only to decrease at the end of 
reionization when the cosmic mean x n approaches zero. The 
anisotropy is most easily observable during the late stages 
where the signal is above 1 mK and is highly anisotropic. 
Note that the anisotropy is stronger in the case of Pop III 
stars, since in that case reionization occurs at a higher red- 
shift and is caused by rarer halos, whose number density 
changes more rapidly with redshift during reionization. 


4 CONCLUSIONS 


At the later stages of reionization, ionization fluctuations 
dominate the 21cm power spectrum (Figures [2] and [3]) and 
the light travel-time delay generates a strong anisotropy. The 
signal is around a few mK on scales and redshifts where 
the time-delay anisotropy is large (Figure [4J, while noise 
and foreground-subtraction should still allow even first- 
generation 21-c m experiments to reach a sensitivity ~ 1 mK 
on large scales feowman et al.l [20061 : iMcQuinn et alJ l2006). 

The strong line-of-sight anisotropy in the ionization 
fluctuations arises in our model since the cosmic mean neu¬ 
tral fraction changes rapidly with redshift and its fluctua¬ 
tions are highly non-linear. The scale and angular depen¬ 
dence of the light-cone anisotropy must be modelled care¬ 
fully in order to separate t he inflationary initial co nditions 
from astrophysical effects llBarkana fe Loe9 2005a ). Large - 
scale numerical simulations of reionization ( Iliev et ahll2006l ) 
could provide guidance for this difficult task, but fully self- 
consistent simulations with hydrodynamics and radiative 


transfer are required. Alternatively, the separation of the 
“physics” (i.e., the inflationary initial conditions) from the 
“astrophysics” may require observations at higher redshifts, 
when the neutral fraction is closer to unity and the 21cm 
power spectrum is dominated by density and peculiar veloc¬ 
ity fluctuations. 
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